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Abstract 
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^ I Diffraction of multi-level atoms by an evanescent wave reflective diffraction grating 

' is modeled by numerically solving the time-dependent Schrodinger equation. We 

. are able to explain the diffraction observed in experiments with metastable Neon. 

I This is not possible using a two-level atom model. The multi-level model predicts 

\^ ' sensitive dependence of diffraction on the laser polarization and on the intensity 

0^ ' ratio of incoming and reflected laser beams. 
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One of the goals of atom optics is the creation of efficient analogues of opti- 
cal elements for atoms. Considerable progress has been made demonstrating 
mirrors, lenses, and beamsplitters [1] . These elements might be combined into 
devices such as atomic interferometers [2-4]. 

Evanescent wave atomic mirrors have been demonstrated in a number of ex- 
periments [5-12]. A reflective diffraction grating for atoms was reported by 
Christ et al. [13]. They observed diffraction of a slowed metastable neon beam 
from an evanescent optical grating formed by counterpropagating laser beams. 
Up to 3% of incident atoms were diffracted by up to 50 mrad from the main 
reflected beam. This constitutes a large angle reflective beamsplitter. Diffrac- 
tion of fast metastable neon atoms has recently been observed by Brouri et 
al. [14]. 

Diffraction is possible because grazing incidence of the atomic beam produces 
an atomic de Broglie wavelength perpendicular to the grating which is com- 
parable to the grating periodicity. The diffraction angles are determined by 
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energy and momentum conservation [12.15]. However to determine the fraction 
of atoms which are diffracted the interaction of the atoms with the evanescent 
field must be analyzed in detail. 

Two-level atom models have dominated the theoretical work on evanescent 
wave devices [15-18]. Although two-level models predict diffraction [17] we 
have previously reported [19] that they cannot explain the particular diffrac- 
tion observed by Christ et al. [13]. In this paper we show that multi- level 
atoms can explain the observed diffraction. It can be understood using the 
quasipotential theory of Deutschmann et al. [17]. The richer structure of the 
multi-level theory allows avoided crossings between ground state quasipoten- 
tials that arc impossible in the two-level theory. These produce diffraction by 
allowing atoms to exit the evanescent field in different quasipotential eigen- 
states than they entered the field [20] . 

The multi-level model not only explains the atomic diffraction observed by 
Christ et al. but also suggests further work. In particular we find a sensitive 
dependence on the polarization of the laser beams [21], and on the incoming 
atomic Zeeman level. For pure s-polarisation, the system can be represented 
as a series of non-interacting two-level systems and hence we do not see any 
diffraction. 



2 The model 

Our numerical model incorporates the ten magnetic sublevels of the 3s [3/2] 2 

^ 3p'[3/2]2 transition of metastable neon. By numerically solving the time 
dependent Schrodinger equation we find diffracted beam fractions consistent 
with the observations of Christ et al. [13]. 

The geometry of the diffraction experiment is shown in Fig. 1. The direction 
perpendicular to the quartz surface is y and the direction parallel to the surface 
X. The origin of the y coordinate is at the interface and the positive direction 
is into the vacuum. We assume the atom is initially propagating towards the 
interface in the xy plane with a positive x-component of velocity. 

The modulus of the evanescent field wavevector Q and its inverse decay length 
q are [24] 



where k is the free field wavevector modulus, assumed to be the same for each 
laser, riq is the refractive index of the quartz, and 6 is the angle of incidence 



Q — kuq sin ^, 




(1) 
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Fig. 1. The geometry of the lasers and atomic beam in the evanescent field diffraction 
experiment. The atomic beam angles are exaggerated for clarity. 

Table 1 

Parameters used in experiment of Christ et al. [13]. 



Parameter 


Value 


Atom 


Ne* 


Atom mass 


3.3 X 10-26 kg 


Atomic velocity 


25 ms-^ 


Angle of incidence 


36 mrad 


Laser detunings 


900 MHz 


Transition wavelength 


594.5 nm 


k 


1.058 X 10^ m-^ 


Q 


1.10 X 10^ m"^ 


q 


2.72 X 10^ m-i 



of the lasers with the quartz surface. The values of these parameters for the 
experiment we shall model are given in Table 1. 

We assume that the atoms arc initially in an cigcnstate \ kox) of the x-componcnt 
of their momentum with eigenvalue hkox- Then because the photon momentum 
of the copropagating (counterpropagating) field is hQ {—hQ) the a;-componcnt 
of the atomic momentum is restricted to the eigenvalues h{kox + nQ), with 
n any integer. We denote the corresponding set of allowed centre-of-mass x- 
momentum eigenstates by {|n),n G integers}. A ground (excited) state atom 
has n even (odd). 

The 3s [3/2] 2 ^ 3p'[3/2]2 transition of metastable neon is between states with 
total angular momentum J=2. Hence there are five ground \mg), and five 
excited |me), sublevels labeled by the magnetic quantum numbers m^/e G 
{—2, —1, 0, 1, 2}, Fig. 2. The discrete basis states for the x-component of mo- 
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Fig. 2. Schematic diagram of the atomic level structure used in the model. The 
squares of the Clebsch-Gordon coefficients C'm^^j^^ for the transitions are shown. 

mentum and the internal state of the atom are therefore 

\n, m) = \n) ® \m). (2) 



Given our initial condition of ground state atoms the magnetic quantum num- 
ber m is rUg (me) if n is even (odd). The atomic motion in the y-direction, 
perpendicular to the quartz surface, is represented by one coordinate basis 
wavefunction '^n,m{y^ t) for each discrete basis state \n, m). We use the follow- 
ing ansatz for the complete state of the atom 

\^{t))^eM-iklj/2Mh)J2'^n,m{y^t)\^^^)^ (3) 

m,n 



where M is the atomic mass. The exponential prefactor accounts explicitly for 
the initial atomic kinetic energy in the x-direction. 

The Hamiltonian H for the system can be divided into two parts: the y- 
component of the kinetic energy Py/2M, with continuous eigenvalues, and the 
rest {V + p'^/2M}, which has discrete eigenvalues, 



Py (Px) is the ^/-component (x-component) of the atomic momentum perpen- 
dicular (parallel) to the quartz surface. V is the sum of the atom's internal 
energy and the electric dipole interaction energy between the evanescent field 
and the atomic transition Heb, 

V= hAa\m,n){m,n\ + Hed, (5) 

m,n(odd) 



This Hamiltonian is in an interaction picture with the atomic dipoles rotating 
at the average of the laser frequencies cu — {cui + u!2)/2, where a;i(ci;2) is the 
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frequency of the laser copropagating (counterpropagating) with the atoms. 
These frequencies can be different, as in the experiment of Stenlake et al. 
[12]. In that case they must be sufficiently similar that the approximation of 
equal photon momentum magnitudes holds. The atomic detuning is the 
difference between the (degenerate) atomic transition frequencies and uj. The 
dipole interaction energy is 

i^ED = d+ • E+ + H.c. (6) 

where H.c. means Hermitean conjugate and d+ is the positive frequency part 
of the transition electric dipole moment operator d 

J 

d+ = ^ \me){me\d\mg){mg\. (7) 

mg,me=—J 

E+ is the positive frequency part of the total electric field 

E+ = SiEiiy, t) e^^iiQx) + £2-^2(2/, t) exp{-iQx), (8) 
Ei{y, t) = exp(-gy) ex.p{-iAdt)Eio. (9) 

Eio {E20) is the amplitude of the copropagating (counterpropagating) field 

at the quartz surface, and is assumed real for convenience. The £j are the 
corresponding field polarization vectors. = {UJ1 — UJ2) /2 is half the frequency 
difference between the two lasers. In the examples we consider in the next 
section A^ = 0. 

We take the z direction as the quantization axis and work with a spherical 
basis of polarization vectors {u±,uo} defined in terms of the Cartesian unit 
vectors by 

U± = ^(Ux ± i^y), Uo = U^. (10) 

The Wigner-Eckart theorem [25] allows us to express the dipole matrix ele- 
ments in terms of Clebsch-Gordon coefficients ^ and a reduced matrix 
element V 

(me|d-u,K) = C^^,„P, (11) 

where p is 0, ±1 for tt, a± transitions. The values of the non-zero Clebsch- 
Gordon coefficients are indicated in Fig. 2. The dipole interaction Hamiltonian 
Eq. (6) then becomes 
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p,mg,me n{even) 



n\ 



+E2{y,t) Y CPn.g,m.e2,p\m„n+l){me,n\ \ +R.C. (12) 

n{odd) ) 



The £i/2,p = • Up are the spherical polarization components of the two 
lasers. 

Substituting the ansatz Eq. (3) into the Schrodinger equation corresponding 
to the Hamiltonian Eq. (4) gives the set of coupled Schrodinger equations 



-1: r 



2Mh 

+ {'r>','me/g\V\n',m'^/g)<i/n',m'^^^ , 



e/g 



h 



S„ = —{2ko,.nQ + n'Q'). 



(13) 

(14) 



We have solved this set of partial differential equations numerically by the 
well known split operator method [19]. Our numerical solutions were checked 
for accuracy by demanding that changing either the time-step or the spatial 
grid size did not significantly alter the final result. 



3 Diffraction 



Before presenting numerical solutions of the Schrodinger equations we consider 
the quasipotentials of the Hamiltonian Eq. (4) [17]. These are the eigenvalues 
of that part of the Hamiltonian having discrete eigenvalues, namely {V + 
pI/2M} Eq. (5). Since V includes the interaction with the evanescent field 
the quasipotentials are a function of y, the distance from the quartz surface. 
The eigenstatcs corresponding to the quasipotentials may be thought of as 
atomic states doubly dressed by the two evanescent fields. 

The eigenstates of V corresponding to the quasipotentials are adiabatically 
followed by atoms moving sufficiently slowly towards the surface, except near 
avoided crossings of the quasipotentials. Hence they behave like actual po- 
tentials, for example slowing atoms down as they climb them, which is the 
origin of atomic reflection [26] . At the avoided crossings non-adiabatic transi- 
tions between the quasipotentials may occur. These make diffraction possible. 
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since the atoms may leave the evanescent field in a different superposition 
of quasipotential eigenstates than that in which they entered, and different 
eigenstates may have different momenta. 

A two-level atom model allows diffraction for certain ranges of parameters. In 
particular the diffraction considered by Deutschmann et al. [17] occurred for 
much lower atomic detunings than were used in the experiment of Christ et al. 
[13]. The lower detunings would have produced unacceptable levels of atomic 
excitation and hence of spontaneous emission. 

Since excited states spontaneously emit, diffraction into ground states, with 
even n, is of most interest. The optimum probability for the n = ±2 diffraction 
orders was estimated to be about 6% by Deutschmann et al. [17]. This is 
because at least four avoided crossings are required, with an optimal transfer 
probability of 50% at each, and 0.5"^ ~ 0.06. The quasipotentials for a two- level 
model of the experiment are shown in Fig. 3(a). Note that the incoming n = 
quasipotential only has avoided crossings with high order quasipotentials, the 
first being with the n — 21 quasipotential [20]. The transition between the 
corresponding eigenstates would involve a 21 photon process and hence be 
very weak, giving negligible n = —2 order diffraction. This is inconsistent 
with the observations of Christ et al. [13]. 

The quasipotentials for the multi-level atom model of the experiment are 
shown in Fig. 3(b). The multi- level model makes possible Raman transitions 
between the ground Zeeman levels. Hence there is now an avoided crossing 
between n = and n = —2 quasipotentials, so diffraction into an n = —2 
eigenstate is expected. 

The experiment of Christ et al. [13] that we have modelled reported up to 
3% diffraction. The ratio of the intensities of the copropagating to counter- 
propagating waves was fixed at 1.64. To calculate the diffracted fraction we 
numerically solved the Schrodinger equations (13), as described in the previ- 
ous section. Our results are quite sensitive to the polarizations of the two laser 
beams [21]. Fig. 4 shows the percentage of diffraction into the n = —2 order 
as a function of the polarization angle of the lasers. The results in this figure 
were calculated for an equal mixture of magnetic sublevels in the incoming 
atomic beam. 

We did not calculate the diffraction from all initial rrig states for all polarization 
angles. This was because each point on Fig. 4 took approximately 30 minutes 
to compute on a VPP300 supercomputer. Typically the spatial grid in the y 
direction had 2048 elements and 13 n states, n = — 6 to n = 6 were used. The 
integration over 1.6 was performed with 0.2 ns timesteps. These numerical 
parameters were varied to ensure the insensitivity of the solutions. 

However we did calculate for all nig states for several representative polari- 
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Fig. 3. Quasipotentials versus perpendicular distance from the quartz surface, for the 
parameters of Table 1. (a) Two-level quasipotentials. (b) Multi-level quasipotentials 
for a laser polarization angle of 5° away from p-polarization. 

sation angles, and in all cases it was found that they produced amounts of 
diffraction bearing a fixed ratio to each other. This is because only one of the 
n = quasipotentials is involved in the diffraction process, and the various 
rUg states enter this quasipotential in fixed ratios. This observation allowed us 
to deduce the overall diffraction for a mixture of magnetic sublevels. 

We have also modelled the conditions which gave maximum diffraction accord- 
ing to unpublished experimental data [20]: The polarization of the copropa- 
gating beam rotated by 5 degrees from perfect p-polarization (electric field 
in the plane of incidence) and the counterpropagating beam rotated 15 de- 
grees from perfect p-polarization [20]. Despite uncertainty as to the accuracy 
of these figures, we found that these parameters gave a large degree of diffrac- 
tion (about 14%). Furthermore, a ratio of 1.64 between the co-propagating and 
counterpropagating beams was cited as an important condition for diffraction 
to occur [13]. Computationally, we found that this condition produced a local 
maximum in the amount of diffraction. 

These results demonstrate that a multi-level model is able to account for the 
diffraction observed by Christ et al. [13]. The difference between the diffraction 
we find (14% for experimental parameters) and that observed, 3%, is reason- 
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Probability for n = -2 diffraction vs. laser polarization 

20 1 , 1 , 1 




Angle of linear laser polarization (°) : 0°= p-polarization 

Fig. 4. Percentage n = —2 order diffraction versus polarization angle for an incoming 
iTig = 2 atom. The percentage of diffraction is plotted versus the linear polarization 
angle, in degrees, of the two laser beams from pure p-polarization, 0°. The solid 
curve is a fit to the calculated points. Other parameters are as in Table 1. 

able given the non-ideal aspects of the experiment and the lack of delicate 
control over the parameters. 

One of the advantages of numerically solving the time-dependent Schrodinger 
equation is that "movies" of the wavefunction evolution are available. The 
wavefunction can be visualized in either coordinate or momentum space. Fig. 
5 shows an initial rUg = 2,n — probability density evolving to produce 
n = —2 order diffraction. 



4 Conclusion 

We have quantitatively modeled, from first principles, the reflection grating 
atomic diffraction experiment of Christ et al. [13]. Our results are consistent 
with the levels of diffraction observed in that experiment. The model predicts 
strong dependence on the polarization of the copropagating and counterprop- 
agating laser beams. This is physically reasonable since the atomic Zeeman 
structure makes the shifts of the various atomic transitions polarization de- 
pendent. There is some evidence for this effect in the experiment [20]. 
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(a) position space 



(c) momentum space 
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(b) position space 
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Fig. 5. Probability densities in position and momentum spaces for an initial nig = 2 
state atom. The probability density Pn = J2m l^n,rraP is plotted against the diffrac- 
tion order n and either the distance from the glass in microns (a) and (b), or the 
velocity perpendicular to the glass (c) and (d). In each case the top figure is the 
initial condition and the bottom figure is 1.6 fis later, after reflection. There is 
no probability in the orders not shown. Parameters are as in Table 1 with laser 
polarizations of 5° away from p-polarization. 



Computational modelling of experiments, such as we have reported, is partic- 
ularly useful if it can suggest new experiments. Our results show that control 
of the polarization of both laser beams is crucial in reflection grating atomic 
diffraction experiments. Furthermore we found that under the conditions of 
the experiment of Christ et al. [13] the Zeeman nig = 2 state produced most 
of the n — —2 order diffraction. Hence optical pumping into this state could 
potentially increase the diffraction. 

Our model is easily adapted to other atoms and we plan to use it to model 
a Caesium evanescent wave diffraction experiment currently underway in our 
group. 
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